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Abstract. The trapped magnetic field is examined in bulk high-temperature 
superconductors that are artificially drilled along their c-axis. The influence of the 
hole pattern on the magnetization is studied and compared by means of numerical 
Q ' models and Hall probe mapping techniques. To this aim, we consider two bulk YBCO 

O . samples with a rectangular cross-section that are drilled each by six holes arranged 

either on a rectangular lattice (sample I) or on a centered rectangular lattice (sample 

11). For the numerical analysis, three different models are considered for calculating the 

trapped flux: (i), a two-dimensional (2D) Bean model neglecting demagnetizing effects 

^^ I and flux creep, (ii), a 2D finite-element model neglecting demagnetizing effects but 

OO ■ incorporating magnetic relaxation in the form of an _E — J power law, and, (iii) , a 3D 

finite element analysis that takes into account both the finite height of the sample and 
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'/^ ' flux creep effects. For the experimental analysis, the trapped magnetic flux density is 
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measured above the sample surface by Hall probe mapping performed before and after 
the drilling process. The maximum trapped flux density in the drilled samples is found 
to be smaller than that in the plain samples. The smallest magnetization drop is found 
for sample II, with the centered rectangular lattice. This result is confirmed by the 
numerical models. In each sample, the relative drops that are calculated independently 
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j^ ' with the three different models are in good agreement. As observed experimentally, 

the magnetization drop calculated in the sample H is the smallest one and its relative 

value is comparable to the measured one. By contrast, the measured magnetization 

drop in sample (1) is much larger than that predicted by the simulations, most likely 

because of a change of the microstructure during the drilling process. 



1. Introduction 

Single crystals of YBa2Cu307_5 containing a series of holes that are artificially drilled 
along a given crystallographic direction have been recently introduced in order to 
improve the performances of YBCO trapped field magnets [H [21 El IH |5] . For instance, 



it was demonstrated that the maximuin trapped flux density in YBCO magnets may 
be enhanced in the presence of holes that favor the oxygen annealing and thus yield 
larger critical current densities [HI [71 El E]. Moreover, when subjected to variable 
magnetic fields, the drilled trapped field magnets exhibit a better thermal stability 
than non-drilled magnets, since they offer a larger surface for exchanging heat with the 
environment. As a result, losses associated with the thermal motion of vortices are 
reduced [10]. Last, by impregnating the holes with resin, the drilled magnets can be 
reinforced against the strong Lorentz forces induced in the bulk of the material [H] . 

As regards the magnetic properties, the holes of a drilled sample influence the 
current stream lines — which cannot cross them — and thus affect the magnetic 
field distribution. In previous works, we modelled the current stream lines in the 
presence of holes and calculated the trapped field in drilled samples of finite or infinite 
heights [12l [13]. It was shown that, for a given critical current density, the maximum 
of the trapped magnetic field is always smaller in a drilled sample than in a non-drilled 
one. Moreover, it was demonstrated that for a given sample geometry, there is an 
optimal position of the holes which minimizes the drop of trapped flux. The holes must 
be arranged in such a manner that their centre lies on the discontinuity line^ of the 
neighboring holes. 

In the literature, although it is has been already demonstrated that drilling holes in 
a sample decreases its trapped flux [TH [15] , no experimental study has been reported so 
far about the influence of the hole pattern on the trapping properties of drilled samples. 
The purpose of this paper is to show experimentally that, considering a given number of 
identical holes, their geometrical arrangement has an influence on the maximum trapped 
flux density. These measurements are confronted to the theoretical predictions made in 
Refs. [121113]. 

The paper is organized as follows. Section [2] describes the preparation and the 
properties of the samples under study. In Section [3l we describe the numerical 
models of increasing complexity (taking progressively into account flux creep and 
geometrical effects) to predict the trapped flux in the samples. Section [4] reports on 
the measurements of the trapped flux before and after the drilling of the holes, and 
Section [5] discuss them in the light of results obtained with the numerical models. 
Section [6] presents the conclusions. 

2. Sample preparation 

In most cases considered so far [3 [3, [HI [9] , drilled samples were obtained from Y211- 
preform where the holes had been arranged before the synthesis of the Y123 phase. This 
approach makes the comparison of drilled samples with different hole arrangements 
difficult, since the critical current densities may differ appreciably from one sample 
to another even though they are processed similarly. We thus followed a different 

I Discontinuity lines are lines near a given hole and across which the current density changes abruptly 
its direction. 
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Figure 1. (1)- Extraction of two rectangular YBCO samples out of a cylindrical one, 
with the help of a wire saw. (2)- Sketch of the hole patterns drilled in the samples and 
pictures of the samples after the drilling process. 



approach and considered bulk YBCO samples that were drilled after their synthesis. 
This approach allowed us to compare the trapped magnetic flux in every sample before 
and after drilling the holes. Provided the microstructure of the sample has not been 
affected by the drilling process, one can make the reasonable assumption that the critical 
current density is unchanged and consider the sole effect of the modification of the 
current stream lines on the magnetic properties of the sample. 

The plain samples were extracted from a top-seeded melt-grown YBCO bulk 
cylinder that was synthesized at CRISMAT (Caen, France). The cylinder has a diameter 
of 22 mm and a height of 16 mm. With the help of a wire saw, two rectangular samples 
were cut out as sketched in Figure [TJ The resulting samples have a length of 8.6 mm, a 
width of 6.6 mm, and a height of 4 mm. 

In each sample, six holes with diameter 1 mm were drilled along the c-axis. To 
prevent the superconducting ceramics from breaking during the drilling process, an 
ultrasonic drilling machine was used. In order to avoid effects associated with the finite 
side lengths of the samples, the six holes were made in the trapezoidal sector delimited 
by the gray dashed lines shown on Figure [H For each sample, this sector lies in the 
rectangular cross section and is dehmited by one long edge of the cross section, a segment 
of its longest median, and two lines which start at the corners with a 45 degree angle. 
This choice is motivated by the critical state configuration found in a non drilled sample 
subjected to a magnetic field applied along the c-axis. In this case, the two boundaries 
at 45 degrees are places where the current stream lines change abruptly their direction, 
whereas they flow parallel to the long edge within the trapezoidal sector [16]. Thus, the 
holes are drilled in a sector where the current lines are initially parallel to one another, 
a situation which is close to that considered in the previous theoretical models [12] . 

For each sample, two series of holes (containing 4 and 2 holes, respectively) are 
made along parallel lines. The hole separation along a line is 1.6 mm and the lines are 
separated by 1.4 mm. The two holes on the second line can be arranged in two different 
ways: (I) either they are aligned with the holes of the first line and form the beginning 
of a rectangular lattice, or (II), the holes are shifted with respect to those of the first 



line by half a separation (0.8 mm) and start forming a centered rectangular lattice. In 
this latter case, the centres of the holes of the second line are in fact placed on the 
discontinuity lines associated with the holes of the first line |12]. A picture of the two 
samples obtained after the drilling process is shown in Figure [H 

3. Numerical models 

The trapped magnetic flux density is calculated in samples I and II by means of three 
different numerical models of increasing complexity: 

(1) A critical state model — or Bean model — 

In the first model, which is a generalisation of the Bean model for arbitrary cross 
sections pTf [18] , the samples are assumed to have an infinite height and a uniform 
critical current density, Jc- The gradient of magnetic fiux density satisfies the 
relation 

(iB J ±/io Jc in penetrated regions, , . 

di \ in virgin regions, 

where P is a point at a given location, i is the total length crossed by the external 
flux to reach P, and fio is the magnetic permeability. This relation can be integrated 
numerically under the constraint that the penetration length i be minimal at each 
point P (or in other words, the shortest penetration path be selected for each 
point P) [12]. The magnetization can then be deduced from the resulting field 
distribution. 

(2) 2D finite element model 

The second model also assumes that the samples have an infinite height but includes 
flux creep effects via the constitutive E — J law 

where Ec is the critical electric field and n is the critical exponent. The field 
distribution can then be solved by integrating the Maxwell equations. The model 
is solved numerically by a finite-element method (FEM) implemented in the open 
source solver GetDP [19]. The detailed procedure has been described in Ref. [13] 
and, for large values of n, uses the properties of a slow magnetic diffusion to reduce 
the number of time steps. In particular, the trapped flux can be calculated with 
only two time-steps: during the flrst step, the applied magnetic flux density is 
increased with a constant sweep rate to a maximum value, it then decreases to zero 
with the same sweep rate during the second step. 

(3) 3D finite element model 

The last model is based on the same equations as those of model (2), but solves 
them for a three-dimensional sample with a flnite height. Details of the model and 
its numerical computation are found in Ref. [13]; in particular, the two time-steps 
method is also used. 



Models (1) and (2) are 2D models that ignore the finite height of the samples. 
Such models are efficient in reproducing the trapped magnetic flux in the median plane 
of a sample, provided its height is larger than the characteristic length of its cross- 
section [201 HI]- This condition is not quite fulfllled for the samples under study. 
Moreover, the 2D models only simulate the trapped flux in the median plane, while 
we only have access to the value at the surface with experiments. However, these 
models are believed to provide a flrst approximation of the relative magnetization drop 
in drilled samples in a much faster way than the 3D model does, without requiring 
excessive calculation loads. These are the reasons why the 2D models are considered in 
this paper. 



4. Measurement of the trapped magnetic flux density 

The trapped magnetic flux density above the sample surface was measured in a Hall 
probe mapping experiment, by moving a miniature probe flxed to a motor-driven xy 
micro-positioning stage over the sample surface [6]. The Hall probe is sensitive to the 
component of the local fleld which is perpendicular to the surface. Its active area is 
0.05 X 0.05 mm^ and its distance from the sample surface is 0.5 mm. The Hall probe 
was moved across the top surface with a step size of 0.5 mm in x and y directions. 
In addition, we also measured the maximum trapped magnetic flux density directly on 
the surface. We brought the Hall probe in close contact with the sample and moved it 
slightly around the sample centre until the maximum reading was found. 

The samples were magnetized with a fleld-cooling process during which the samples. 
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Figure 2. Distribution of the trapped magnetic flux density at 
sample surface of samples I and II, before and after drilling. 



0.5 mm above the 
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initially at room temperature, were immersed in liquid nitrogen in a uniform magnetic 
flux density of 300 mT created by a large copper coil. The field-cooling process lasted 
5 min, and the characterization of the trapped field started 15 min after the applied 
magnetic field had been switched off, so that magnetic relaxation effects were negligible. 
The distribution of the trapped magnetic flux density at 0.5 mm above the top 
surface of sample I — left panels — and sample II — right panels — are shown in 
Figure [21 respectively before drilling, (a), and after drilling, (b). In sample I, the trapped 
magnetic flux density has a single maximum at -B™appcd ~ 40 mT before drilling; after 
drilling, the maximum field drops to -Snapped ~ 25 mT. In sample II a larger maximum 
trapped flux density is found before drilling, with -B™appcd ~ 130 mT; it reduces to 
-^trapped ~ 100 '^T after drilling. The relative drop of the maximum trapped flux 
density is found to be larger in sample I. It can also be observed that, in addition to 
reducing the trapped flux, holes also modify its distribution with the largest effects near 
the maximum. 

5. Discussion 

The measured drops in trapped flux of the drilled samples are compared to the simulation 
results. The same set of parameters are used for drilled and plain samples. They are fixed 
as follows. First, the critical current density is chosen such that model (3) reproduces 
well the value of the maximum trapped flux density that is measured on the top surface 
of each plain sample. This gives Jc = 4.1 10^ A/m for sample I and Jc = 8.8 10^ A/m 
for sample II. Second, for the E — J power law in models (2) and (3), the critical 
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Figure 3. Two-dimensional trapped magnetic flux distribution at « 0.5 mm above 
the sample surface of the drilled samples I and II. The top panels correspond to 
the simulation results with the 3D-FEM models and the bottom panels to the 
measurements . 



Table 1. Comparison of the maximuni trapped magnetic flux density in samples I 
and II 



Sample I 


Median plane 
Bean model FEM - 2D 


FEM - 3D 


Top 
FEM - 3D 


SURFACE 

Measurement 


Before drilling 
After drilling 


169 mT 
126 mT 




137 mT 
104 mT 


95 mT 
75 mT 


61 mT 
46 mT 


60 mT 
33.7 mT 


Relative drop 


25% 




24% 


21% 


25% 


44% 




Sample II 


Median plane 
Bean model FEM - 2D 


FEM - 3D 


Top 
FEM - 3D 


SURFACE 

Measurement 


Before drilling 
After drilling 


358 mT 
291 mT 




310 mT 
253 mT 


244 mT 
207 mT 


154 mT 
130 mT 


155 mT 
120.7 mT 


Relative drop 


19% 




18% 


15% 


16% 


22% 



exponent n is taken a,s n k, 25, a typical value for YBCO bulk samples [22] . 

Unlike the experiments where the samples are magnetized with a field-cooling 
process, the simulations reproduce a zero field-cooling process, and hence require a 
maximum applied field larger than 300 mT as used in experiments. In order to fully 
magnetize the samples, even those with an infinite height, the maximum value of the 
applied magnetic flux density is chosen at 1 T and the sweep rate is equal to 1 mT/s. 
Let us first compare the distributions of the trapped flux above the surface of the drilled 
samples with that calculated with the 3D FEM model. Figure E] shows the measured 
trapped flux distribution at 0.5 mm above the surface in the bottom panels and the 
numerical predictions in the top panels, for each sample. Note that in the experimental 
results, the holes cannot be resolved because of the distance between the Hall probe and 
the surface of the samples. In sample I, the simulated distribution of the trapped flux 
is symmetric, with a maximum trapped flux density occuring at the center of the cross 
section. This distribution deviates from the experimental mapping (bottom left panel), 
where the maximum trapped field is found to be slightly shifted to the left of the centre. 
The agreement between the calculated and the measured distributions is much better 
for sample II, where the maximum of the trapped flux is found to be shifted to the right 
of the centre in both the experimental mapping and the simulated distribution. 

The measurement of the trapped flux distribution at a finite distance above the 
sample surface can be affected by a systematic error in the distance between the Hall 
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Figure 4. Micrographs of tire bottom and top surfaces of sample I and II after the 
drilling of the holes. 



probe and the sample surface that cannot be reproduced accurately in simulations. It 
is thus interesting to compare the maximum trapped flux densities when the probe is 
in contact with the sample. Tabled] reports the simulated and the measured maximum 
trapped flux density in sample I and II, before and after drilling. 

The maximum trapped flux density is simulated in the median plane of the samples 
with models (1), (2) and (3) and on the surface with model (3). It is found that the 
maximum trapped flux density is in each case the largest with the Bean model, the 
flux creep and geometrical effects bringing additional drops of trapped flux. Further, 
although the drops in magnetic flux density are of different magnitudes in each model, 
the relative drops are found to be nearly equal, about 23 % for sample I and 17% for 
sample II. These results indicate that the centered rectangular lattice is the conflguration 
with the smallest drop in the trapped flux, as was already found in samples with different 
cross sections in Ref. |12] . 

Experimentally, the measured maximum magnetic flux density is found to drop 
by 44% in sample I, and by 22% in sample II. On the one hand, the smallest relative 
drop that is observed in sample II is well reproduced by simulations, which conflrm the 
previous theoretical results in Refs. fi2[ [T3] . It is interesting to note that the simple 
2D models, while not absolutely suitable for samples under consideration in this study, 
already yield a reasonable estimate of the relative drop that is measured at the surface of 
sample II. For sample I, on the other hand, the numerical simulations fail in describing 
the large measured drop. 

Given the discrepancy between theory and experiment for sample I, the large 
observed fleld drop cannot be only due to a reorganisation of the current lines after 
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drilling. Presumably, drilling also affected the microstructure of the sample. Figure H] 
shows optical microscope pictures of top and bottom surfaces of sample I and II after 
drilling. The arrows (a) and (b) show a crack between two neighboring holes in sample 
I. This crack is found on both surfaces and has probably be produced during drilling. 
Its position is consistent with the shift to the left observed in the maximum of the field 
distribution in Figure |3l By contrast, sample II presents a crack on one of its sides 
(arrow (c)) but not on the other, so that drilling appears to have preserved most the 
microstructure of this sample. 

6. Conclusion 

We have shown with experiments and modelling that the arrangement of the holes in a 
drilled sample influences the trapped magnetic flux. We considered two samples with the 
same dimensions and made six holes in each, according to two different arrangements. 
In sample I, the holes are placed on a rectangular lattice, and in sample II, they form a 
centered rectangular lattice. 

The measurements of the trapped flux density above the sample surfaces have shown 
that the drilled samples trap a smaller magnetic flux than a sample with no hole. The 
sample II with the centered rectangular lattice exhibits the smallest relative drop of the 
maximum trapped magnetic flux. 

The measurements have been compared to the results obtained with numerical 
simulations. We have considered three models: (1), a 2D Bean model neglecting the 
flux creep and assuming that the samples have an infinite height, (2), a 2D FEM model 
with an E — J power law, and (3), a 3D FEM model taking into account both the flux 
creep and the finite size of the samples. A good agreement between the simulations with 
each of the three models and the measurements has been found for the relative drop of 
trapped flux in sample II. Sample I exhibits a larger drop than what is predicted by the 
simulations, most probably because drilling produced cracks over the full sample height. 
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